Entropic effects on the Size Evolution of Cluster Structure 
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We show that the vibrational entropy can play a crucial role in determining the equilibrium 
structure of clusters by constructing structural phase diagrams showing how the structure de- 
pends upon both size and temperature. These phase diagrams are obtained for example rare 
gas and metal clusters. 

PACS numbers: 61.46.+w,36.40.Mr,36.40.Ei 



Much of the interest in clusters or nanoparticles de- 
rives from the insights they can provide into how proper- 
ties emerge and evolve on going between the atomic and 
molecular and bulk limits. Cluster structure provides a 
particular interesting example of this size evolution. At 
large enough sizes the clusters must display the bulk crys- 
talline structure, but this limit may sometimes only be 
achieved at very large sizes (e.g. at least 20 000 atoms 
for sodium clusters Q) and before that limit is reached 
unusual structural forms are often observed. For exam- 
ple, many clusters bound by van der Waals or metallic 
forces exhibit structures with five-fold axes of symmetry, 
a possibility that is forbidden in bulk crystalline mate- 
rials. For these clusters the dominant structural motif 
typically changes from icosahedral (Fig. 0(b)) to decahe- 
dral (Fig. 0(c)) to face-centred-cubic (fee) (Fig. 0(a)) as 
the size increases. 

For many materials these structural changes occur at 
sizes that are too large for global optimization to be 
feasible. Therefore, the typical theoretical approach to 
systematically investigating the size evolution of cluster 
structure is to compare the energies of stable sequences of 
structures, such as the forms shown in Fig. 0. 'Crossover 
sizes' are then identified where the sequence with lowest 
energy changes. At this crossover the most common equi- 
librium structure is expected to change. This technique 
has been applied to rare gas g [|, [|, metal || |[ |7[ ||] 
and molecular clusters ||. 

The above approach is certainly valid at zero temper- 
ature, since the equilibrium structure then corresponds 
to the one with lowest energy. At other temperatures, 
however, the structure with lowest free energy needs 
to be found. However, perhaps through an expecta- 
tion that entropic effects are unlikely to be important or 
are too complicated to take into account, size is usually 
the only variable that is considered both experimentally 
@ and theoretically §, § g | g % §. 

In this paper we consider the role that entropy plays in 
the size evolution of cluster structure, and show that tem- 
perature can be a key variable in determining the equi- 



FIG. 1: Three examples of the structures clusters can 
adopt: (a) a fee 38- atom truncated octahedron, (b) a 
55-atom Mackay icosahedron Jl3| ], and (c) a 75-atom 
Marks decahedron G3] . These clusters have the optimal 
shape for the three main types of regular packing seen in 
clusters: face-centred cubic, icosahedral and decahedral, 
respectively. The latter two structural types cannot be 
extended to bulk because of the five-fold axes of sym- 
metry. 



librium structure of a cluster. A clue to this result can be 
garnered from the growing number of examples of solid- 
solid transitions in clusters where the structure changes 
from fee or decahedral to icosahedral as the temperature 
increases |lj| [li], [l?], 18, IjJ. The most well-investigated 
examples of these transitions are for those small Lennard- 
Jones (LJ) clusters that have a non-icosahedral global 
minimum. For N < 150 there are only 8 examples and 
these occur at sizes where the non-icosahedral morphol- 
ogy has a near optimum shape, whereas the icosahedral 
structures have an incomplete outer layer. The global 
minimum is fee at iV=38, decahedral at N=75-77, 102- 
104 and has an unusual structure called a Leary tetrahe- 
dron at N=98 

Here we analyse these examples further in order to 
identify what is the most significant contribution to the 
difference in entropies between the structural types, par- 
ticularly as the size increases. Direct computation of T ss , 
the temperature of the solid-solid transition, is possible 
by parallel tempering |2l| , but this is computationally de- 
manding because of the large (free) energy barriers be- 
tween the structural types Jl7], [23], and has only been 
done for LJ38 p3fl. 
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TABLE I: Estimates of the solid-solid transition temper- 
ature T ss for those Lennard- Jones clusters with less than 
150 atoms that have non-icosahcdral global minima. V a 
and Vq are the mean frequencies of the lowest-energy 
non-icosahedral and icosahedral minima, respectively, e 
is the equilibrium pair well depth of the LJ potential. 
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The superposition method |24[ 25 provides another 
approach to this calculation and one which is particu- 
larly useful for our present purposes because it allows us 
to analyse the different contributions to the entropy. In 
this approach the partition function is written as a sum 
over all the minima on the potential energy surface, and 
by restricting the sum to a certain subset of the min- 
ima the thermodynamic properties of a particular region 
of configuration space can be obtained. The centre of 
the solid-solid transition occurs when the partition func- 
tions for the two competing structural types are equal, 
i.e. Za = Zb- In the harmonic approximation pq] this 
then gives 
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where (3=1/ kT, Ei is the potential energy of minimum 
i, Vi is the geometric mean vibrational frequency, and 
m=2N\/hi is the number of permutational isomers of 
i, where hi is the order of the point group. T ss values 
calculated using Eq. (Q) are given in Table 0. It is note- 
worthy that there is a general decrease in the values with 
increasing size. 

One contribution to the entropy of a morphology comes 
from the number of low-energy minima of that type. 
For the current examples this term always favours the 
icosahedra because there are many low-energy icosahe- 
dral minima with different arrangements of the atoms in 
the incomplete outer layer (22). Another entropic con- 
tribution can come from the symmetry of the cluster. 
At ^=38, 75 and 98 the global minimum has high sym- 
metry, thus reducing the number of permutational iso- 
mers. For these three sizes this factor again favours the 
icosahedra. Lastly, there is the vibrational entropy. For 
the current examples the icosahedral structures have a 
smaller mean vibrational frequency, which again favours 



the icosahedra. However, unlike the previous two contri- 
butions to the entropy, this term favours the icosahedra 
for a LJ cluster of any size. Furthermore, the absence of 
any solid-solid transitions when the LJ global minimum 
is icosahedral therefore suggests that the vibrational en- 
tropy is crucial. 

We can analyse further the effect of the vibrational en- 
tropy by applying an Einstein approximation, i.e. we as- 
sume that all the minima have the same mean vibrational 
frequency. The resulting values for T ss are also given in 
Table Q. Although the transitions still occur, they do 
so at significantly higher temperature with the error in- 
creasing with size, because 77 is raised to the power 3N—6 
in Eq. ([!]). This increasing dominance of the vibrational 
entropy is the main reason that the actual T ss generally 
decrease with size. These results, therefore, suggest that 
solid-solid transitions, rather than being unusual, should 
be expected for systems where different structural types 
have systematic differences in V at sizes where the mor- 
phology with lower vibrational entropy is the lowest in 
energy. 

To obtain our results for T ss we systematically searched 
the low-energy regions of the potential energy surface for 
these clusters in order to generate the relevant samples 
of low-energy minima. However, this is not a practical 
approach at large sizes, therefore we seek a simpler way 
to estimate T ss . Firstly, if we assume that all the minima 
associated with a morphology have the same energy and 
vibrational frequency, then 
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where AE=Eb — Ea and ua is the total number of min- 
ima (geometric and permutational isomers) associated 
with morphology A. Although in all the specific exam- 
ples in Table | there are far more icosahedral low-energy 
minima, on average the number of minima will be ap- 
proximately the same. Furthermore, at large sizes the 
vibrational term will dominate the denominator. In this 
limit 
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We tested this expression for the examples in Table | us- 
ing the properties of the lowest-energy minima of the two 
competing morphologies. As expected this estimate be- 
comes more accurate as the size increases, and the errors 
for the largest sizes are small. 

We can also use Eq. (^) to determine how crossover 
sizes depend on temperature. As with the theoretical 
studies mentioned earlier, we can compare stable se- 
quences of structures, but now monitoring not only their 
energies but also their vibrational frequencies. These 
properties can be fitted to the forms 



E = a E N + 6 s iV 2/3 + c E N 1/3 + di 
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FIG. 2: (a) Energies of the competing structural types 
for LJ clusters. The data points represent clusters with 
the optimal shape at that size, and the continuous lines 
are fits to these data using Eq. (||). The energy zero is 
£toct) a fit to the fee truncated octahedra with regular 
hexagonal {111} faces (Fig. |l|(a)). (b) Structural phase 
diagram for LJ clusters. The data points represent the 
melting temperatures of the four smallest Mackay icosa- 
hedra. For argon lefc _1 =121K. 
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where the first two terms represent volume and surface 
contributions, respectively. These expressions can then 
be input into Eq. (||) to map out the solid-solid tran- 
sitions in the structural phase diagram. Finally to com- 
plete the diagram, the melting line has to be determined. 
Again it will have the form 
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We first present results for LJ clusters, which provide 
a reasonable model of heavy rare gas clusters, such as ar- 
gon. Up to A«10 000 the energies of the stable sequences 
of structures were obtained by minimization of all degrees 



of freedom. Above this size minimizations were possi- 
ble up to A«35 000 for decahedral and A«80 000 for 
fee structures if the cluster was constrained to maintain 
the correct point group symmetry (D5/1 for decahedral 
and Oh for fee). The energies of the optimal clusters 
for each structural type are shown in Fig. |2|(a). The 
lobed shape of the decahedral and fee lines results from 
changes in the shape of the optimal sequence; as the size 
increases the forms depicted in Fig. |l|(a) and (c) become 
more rounded by the introduction of {110} facets of in- 
creasing size, and the grooves in the Marks decahedra 
become deeper. The energetic crossover sizes are con- 
sistent with previous results Q]: Ai cos ^dcca(? 1 =0)wl450 
and A dcca ^ fcc (T=0)«213 000. 

In order to calculate V the Hessian matrix must be di- 
agonalized, and so the sizes for which V can be calculated 
are limited to A<3500. Above this size we rely upon 
the extrapolation of Eq. To construct the melting 
line a m is assigned the value of the zero pressure bulk 
melting temperature |27 and the other three parame- 
ters in Eq. (||) are fitted using the melting points of the 
first four Mackay icosahedra, which were obtained from 
Monte Carlo simulations. 

The structural phase diagram that results from these 
calculations is shown in Fig. ||(b). When interpreting 
these diagrams, one should remember that the phase 
boundaries divide the plane into regions where the ma- 
jority of clusters have a particular equilibrium structure. 
If the non-monotonic variation of cluster properties, as 
illustrated by the examples in Table |[ were fully taken 
into account, the phase boundaries would be considerably 
rougher. 

The effect of the vibrational entropy can be clearly seen 
from the slopes of the phase boundaries in Fig. |2|(b), e.g. 
A r icos-*doca(T 1 m )«6550 and A dcca ^ fcc (T m )wl0 600 000. 
At higher temperatures icosahedra and Marks decahedra 
remain most stable up to considerably larger sizes This 
is because of the relative values of their vibrational fre- 
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are 2.06% and 0.19% less than 



, respectively. Even though the difference between fee 
and decahedral frequencies is much smaller, the effect on 
the slope of the phase boundary is larger because of the 
larger value of N in the denominator of Eq. (||) . 

The techniques that we have developed are readily ap- 
plicable to other atomic clusters. This is illustrated in 
Fig. H for silver clusters described by the Sutton-Chen 
potential Bq| . The phase diagram has a very simi- 
lar form to the LJ diagram except that the crossovers 
occur at smaller sizes, and again the icosahedra and 
decahedra are substantially stabilized by temperature. 
Mcos^dcca(0)«240 and A doca ^ fcc (0)«36 600 compared to 
-^Vicos^dcca(Tm)«3500 and iVd cca ^f C c (T m )«253 000. 

However, the method cannot be simply applied to 
molecular clusters because the orientational degrees of 
freedom add an extra degree of complexity. Tests on 
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FIG. 3: Structural phase diagram for silver clusters. 



N2 clusters showed that for minima where the molecular 
centres have a common structure but the molecules have 
different orientations there is a considerable variation in 
v. Thus one of the assumptions underlying the use of 
Eq. (||) is invalidated and one would have to consider 
the competition between orientational isomers as well as 
between structural types. 

The key role that we have shown the temperature to 
play in determining the equilibrium structure of a clus- 
ter underlines the importance of only comparing clusters 
that have both the same size and temperature. This may 
help to explain some of the apparent contradictions be- 
tween experiments and zero-temperature theoretical cal- 
culations and between different experiments (e.g. elec- 
tron diffraction @ and EXAFS JO] predict different 
crossover sizes for rare gas clusters). 

Although the developments in this paper allow a more 
complete comparison of theory and experiment, it should 
be remembered that equilibrium may be hard to achieve 
in experiment. For example, experiments on gold clus- 
ters E9] and clusters of Ceo molecules [30 show the im- 



portance of annealing the clusters at sufficiently high 
temperature to locate the most stable forms. Further- 
more, the (free) energy barrriers between structural types 
jl7], |22|, which have been shown to be large for LJ clus- 
ters, are likely to increase with size, making structural 
transformations increasingly hard. In such a situation, 
growth may preserve the dominant structure at the last 
size at which equilibrium was obtained |3l|] , or may occur 
around a core structure that allows rapid growth |32[ . 
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